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ABSTRACT 

Accretion disks around supermassive black holes are widely believed to be the dominant source of 
the optical- ultraviolet continuum in many classes of active galactic nuclei (AGN). We study here the 
relationship between the continuum colors of AGN and the characteristic accretion disk temperature 
(T max ). Based on NLTE models of accretion disks in AGN computed as described by Hubeny et al. 
(2000), we find that continuum intensity ratios for several pairs of wavelengths between 1350 and 
5100 A should show a trend of bluer colors for higher T max , notwithstanding random disk inclinations. 
We compare this theoretical expectation with observed colors of QSOs in the Sloan Digital Sky Survey 
(SDSS), deriving black hole mass and thence T max from the width of the Mg n broad emission line. 
The observed colors generally do not show the expected trend and in some cases show a reverse trend 
of redder colors with increasing T max . The cause of this discrepancy does not appear to be dust 
reddening or galaxy contamination but may relate to the accretion rate, as the offset objects are 
accreting above ~ 30 % of the Eddington limit. The derived disk temperature depends primarily on 
line width, with little or no dependence on luminosity. 
Subject headings: galaxies: active — quasars: general — black hole physics 
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1. INTRODUCTION 

Accretion onto supermassive black holes is widely ac- 
cepted as the energy source of AGN. Much of the optical 
and ultraviolet luminosity, specifically the so-called "Big 
Blue Bump" may be thermal emission from an accre- 
tion disk (Shields 1978; Malkan 1983; for a review of 
attempts to relate disk models to QSO spectra see Ko- 
ratkar & Blaes 1999). In the simplest picture, radiation 
from each radius in the disk is powered by energy pro- 
duced locally by viscous dissipation within the vertical 
thickness of the disk. Models achieve some success in 
fitting the spectral energy distribution (SED) of individ- 
ual QSOs (e.g., Sun & Malkan 1989; Blaes et al. 2001; 
and references therein). Simple disk models nevertheless 
have shortcomings (e.g., Antonucci 1999). Time varia- 
tions of the optical, ultraviolet, and X-ray continuum of 
some AGN show patterns suggestive of reprocessing in 
an externally illuminated atmosphere (Collin 1991). Po- 
larization of the optical continuum of QSOs is weaker 
in degree and opposite in sign from that expected from 
a smooth electron scattering atmosphere (see Coleman 
& Shields 1990, and references therein). The soft X- 
ray continuum cannot be explained by ordinary thermal 
emission from the disk atmosphere, and requires another 
ingredient such as a Compton-scattering corona (Wilkes 
& Elvis 1987; Comastri et al. 1992; Blaes et al. 2001). 
On the other hand, Kishimoto et al. (2005) argue from 
spectropolarimetry that the infrared continuum of AGN 
accretion disks has a slope similar to that expected from 
simple disk models. 

The mixed success of the disk model makes it impor- 
tant to explore further tests of the basic agreement be- 
tween the model and observations. One simple test in- 
volves the temperature of the disk and the expected color 



of the continuum. Wien's law suggests that hotter disks 
should in general have bluer continua. This is born out 
by models, as we discuss below. The disk effective tem- 
perature varies with radius in a self-similar way, and the 
temperature of a given disk can be parameterized by 
the maximum value of the local effective temperature, 

Tmax oc L^Mbh -1 ^ 2 , that occurs at a radius slightly 
outside the inner disk boundary (see review by Novikov 
&; Thorne 1973). T max varies with Mbh and the accre- 
tion rate M as T max 4 oc M/Mbh 2 , because the bolomet- 
ric luminosity ( Lboi) is proportional to Mand the disk 
radius scales with the gravitational radius R g = GM/c 2 . 
One consequence is that, for objects shining at a given 
fraction of the Eddington limit (Lea oc .Mbh)) the tem- 
perature decreases with increasing luminosity or mass. 
On this basis, one might expect luminous QSOs to have 
softer (redder) continuum slopes than lower luminosity 
AGN. This expectation is not consistently fulfilled by 
observed SEDs of AGN. There is a trend of increasing 
soft X-ray excess with narrower H/3 and higher Edding- 
ton ratio (e.g., Laor et al. 1997); these quantities corre- 
late in practice with higher T max . The "Baldwin effect", 
which involves decreasing emission-line equivalent width 
for higher luminosity, is also consistent (Wandel 1999; 
Dietrich et al. 2002). More luminous objects should 
have cooler disks and softer ionizing continua. Scott et 
al. (2004) find harder EUV slopes for higher Mbh and 
luminosity, although curiously the trend is largely lost in 
scatter in their plot of slope versus T max . Blaes (2003), 
looking at continuum slope as a function of Mbh and 
L/L^d for a small set of AGN, found a large degree of 
scatter but an overall trend apparently opposite to the- 
oretical expectation (also noted by Cheng, Shields, & 
Gebhardt 2004). 
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T max may be estimated for individual objects with 
knowledge of AIbh and M. The accretion rate can be 
estimated from the bolometric luminosity, and A/bh can 
be estimated with the aid of recently established methods 
utilizing the widths of the broad emission lines (Kaspi et 
al. 2000, 2005; Bentz et al. 2006). In this paper, we use 
spectra from the SDSS DR3 1 to examine the continuum 
slope of QSOs as a function of T max . In Section [2 we de- 
scribe the theoretical models we use and compare them 
to observed QSOs. Results are discussed in Section [3l 

2. MODELING AND DATA SAMPLING 

2.1. Accretion Disk Properties 

The Eddington limit for a black hole (BH) with mass 
A/ B h is given by 

47tcGM B h 
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where n e is the electron scattering opacity per unit mass 
and M 8 = M BH /(10 8 M ). The disk luminosity produc- 
tion through accretion onto the black hole is 



L = eMc 2 = (10 45 76 erg s _1 )e_iM . 



(2) 



where Mq is the accretion rate in M Q yr _1 , and e = 10e_i 
is the efficiency For a Schwarzschild (non-rotating) black 
hole, e = 0.057, and for a rapidly rotating Kerr hole 
with angular momentum parameter a» = 0.998, e = 0.31 
(Novikov & Thorne 1973). The effective temperature at 
radius r is given by the local radiative flux leaving the 
atmosphere, T cS = {F/a) 1 / 4 . For a* = 0.998 (Thorne 
1974), 
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or alternatively, T max = (10 5 54 K)M 8 1/2 L% A , where 
L 46 = Lboi/(10 46 erg s^ 1 ). For a Schwarzschild hole, 
Tmax is cooler by a factor of 10 046 for a given Mbh and 
Lboi (e.g., Shields 1989). 

In traditional treatments, the local vertical structure is 
determined by hydrostatic equilibrium and radiative and 
convective transport of heat produced locally by viscous 
stresses. The viscosity also determines the radial diffu- 
sion of matter and the surface density of matter in the 
disk as a function of radius. The locally emitted spec- 
trum is determined by T c g and the effective gravity in the 
atmosphere. The locally emitted spectrum shows strong 
features such as the Lyman edges of hydrogen and he- 
lium, but in the observed spectrum these features are 
broadened by relativistic effects. 

2.2. Deriving A/bh and T max 

Derivation of Mbh from the width of the H/3 or Mg n 
broad emission lines has become widely accepted in re- 
cent years. The FWHM of the broad lines is taken to give 
the circular velocity of the broad-line emitting material 
(with some geometrical correction factor). The radius of 
the broad-line region (BLR) , derived from echo mapping 
studies, increases as a function of the continuum lumi- 
nosity, R(x L r with r = 0.5 - 0.7 (Wandel et al. 1999; 
Kaspi et al. 2000, 2005). Simple considerations involv- 
ing photoionization physics suggest T — 0.5 (Shields et 

1 The SDSS website is |http: / /www.sdss.org| 



al. 2003). Bentz et al. (2006) find T = 0.52 ± 0.04 af- 
ter correcting the sample of Kaspi et al. (2005) for host 
galaxy contamination. Such studies typically use the H/3 
broad line, but McLure & Jarvis (2002) find that the 
Mg II A2800 line is a valid alternative useful at higher 
rcdshifts. Here we use the expression for A/bh 

(4) 



M BH = (10 7 - 69 M M OOO L^, 



where v 3000 = FWHM/3000 km s _1 , T = 0.5, 
L 44 = AL A (5 1 00)/(10 44 erg s" 1 ), and the coefficient 
is taken from Shields et al. (2003). We use the FWHM 
of the Mg II line in order to get maximum coverage in 
redshift for a single broad emission line. 

The derived value of T max depends almost entirely 
on the the broad line FWHM for a given object. The 
bolometric luminosity can be estimated as Lboi = /l x 
AL A (5100), following Kaspi et al. (2000). Using Eq. H 
and this expression for Lboi, Eq. [3] becomes 
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Note that for T = 0.5, as used here, the derived disk 
temperature is inversely proportional to FWHM with no 
dependence on L, save the bolometric correction factor 
f L = 9 (Kaspi et al. 2000). Shang et al. (2005), in 
an analysis of the SEDs of 17 well observed QSOs, also 
adopted f L = 9, although their broad band SEDs typi- 
cally indicate a larger value, on average Jl ~ 13. More 
recently, Richards et al. (2006) find an average bolomet- 
ric correction from 5100 A of ~10 from composite SEDs 
of 259 QSOs. We discuss the sensitivity of our results to 
the value of Jl in Sec. 13.11 below. 



2.3. The Color-T„ 



Relation from the Non-LTE 
Models 



One might expect that cooler disks will in general have 
redder colors. This certainly would be true (apart from 
inclination effects) for non-relativistic disks emitting lo- 
cally as black bodies. However, disks in AGN are ex- 
pected to depart from black body emission locally, and 
the energy distribution received by a particular observer 
is strongly affected by relativistic effects. The observer's 
angle i with respect to the disk axis is also important, 
because the radiation from the innermost, hottest part 
of the disk will be relativistically beamed near the equa- 
torial plane, especially for rapidly rotating holes (Cun- 
ningham 1975; Laor & Netzer 1989). Therefore, models 
are needed to establish the expected statistical relation- 
ship between T max and color. In recent years, I. Hubeny 
and his colleagues have computed NLTE disk models for 
AGN and made comparisons with observation (Hubeny 
et al. 2000). Hubeny et al. have kindly made avail- 
able the computer program AGNSPEC which gives the 
observed energy distribution (AitdLv/dVL) emitted in the 
observer's direction for model disks as a function of A/bh, 
Al, a„, inclination, and the Shakura & Sunyaev (1973) 
viscosity parameter a. The vertical structure and local 
emission is evaluated at each radius in the disk, and the 
spectrum viewed at infinity is calculated via the gen- 
eral relativistic transfer function. Using AGNSPEC, we 
computed a set of models with values of Mbh and M 
representative of our SDSS QSO sample (discussed be- 
low). We select the range of Mto give Eddington ratios 
consistent with those in our observed sample (see Fig. [7] 
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Fig. 1— Continuum flux ratio F A (A5100)/F A (A1350) for 
Schwarzschild and Kerr black holes. The Schwarzschild models 
are cooler because they were computed for the same set of values 
of A/bh an d ^bol- The points show an overall trend of bluer colors 
with increasing T ma x- Models with high inclination (cos i < 0.5, 
black dots) show more scatter due to the relativistic effects on 
edge-on disks. 

and discussion below. We fix a = 0.10 and, except as 
noted below, a* = 0.998. Hubeny et al. (2000) show 
that models with a = 0.1 and 0.01 differ negligibly at 
the wavelengths of interest here. To assess the effect of 
disk inclination on the observed colors, we assign a ran- 
dom inclination (uniform distribution in cos i) to each 
model. 

In Fig. Q] we plot 
L A (5100 A)/Z A (1350 A) vs. 
models with cos i > 0.50. 
limit to the inclination angle for observable quasars if 
the "unified" model of AGN, involving an obscuring 
torus, applies to our QSOs (Antonucci & Miller 1985). 
There exists a clear overall trend wherein the spectra 
become more blue with increasing disk temperature, 
consistent with intuition. The models with cos i < 0.50 
add more scatter, since relativistic beaming and Doppler 
shift effects result in more flux at shorter wavelengths. 
However, the overall trend of bluer color with increasing 
Tmax is still evident. We include models for non-rotating 
black holes for the same Mbh an d i-boi as the Kerr 
models. These have lower T max by 0.46 dex for given 
Mbh and Z-boi, but they conform to the same trend of 
color with T max . 
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2.4. The Color-Tmax Relation from the SDSS Datasets 

We examined the color (i.e., continuum flux ratios) 
of SDSS QSOs as a function of T max derived from the 
FWHM of Mg II as in Eq. E (see Sec. E2J. The 
SDSS spectra cover a wavelength range of A3800 to 
A9200 A. The rest-frame colors measured for our sample 
are r 54 = F A (5100)/F A (4000), r 42 = F A (4000)/F A (2200), 
and r 2i = F A (2200)/F A (1350). The objects have red- 
shifts between z ~ 0.38 and z ~ 2.0, such that Mg II lies 
well within the SDSS spectrum. The broad Fe n emis- 
sion blends were removed with the aid of templates in the 
optical and UV bands (Marziani et al. 2003 and Vester- 
gaard & Wilkes 2001, respectively). Continuum fluxes 
were computed by taking the mean (biweight) value of 
the flux over a 30 A window centered on the wavelength 
of interest. The chosen wavelengths are relatively unaf- 
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Fig. 2. — Observed continuum colors for individual QSOs as 
given by flux ratios at three pairs of wavelengths. The abscissa is 
log T max derived as described in the text, on the assumption of a 
rapidly rotating black hole. 



fected by Fe n emission. The FWHM of the broad Mg n 
AA2796, 2803 doublet was measured by means of a least 
squares fit of a Gauss-Hermite function (Pinkney et al. 
2003) to the line profile, together with a linear fit to the 
continuum in the vicinity of the line. Fits were accepted 
if the line widths and equivalent widths had an accuracy 
better than 15% and visual inspections showed a good fit 
and no artifacts in the spectrum or continuum windows. 
For more details on the measurement procedure, see Sal- 
viander et al. (2006). Figure shows the results for indi- 
vidual QSOs in the three wavelength ranges. The figures 
show substantial scatter, but overall a flat or upward 
trend, particularly for higher T max . In order to clarify 
the trends, we have binned the data for each wavelength 
range such that there are approximately equal numbers 
of objects per bin of T max (see Fig.[3J For a given bin, we 
averaged the flux ratios log r^, log r4 2 , and log r%\, for 
each object for which these data exist. At longer wave- 
lengths, agrees fairly well with the predictions of the 
Kerr models. The observed r 42 tracks the models at low 
Tmax but shows a significant upturn above ~ 10 5 4 K. 
Our shortest wavelength pair, 2200/1350 A, disa grccs 
with theoretical predictions even at low disk tempera- 
ture, and the departure increases at higher T max . The 
observed trend of r%i is nearly flat with Tmax, not show- 
ing the 'U'-shape seen in r^.- In order to display a wider 
wavelength range than is possible with SDSS spectra of 
a given redshift, we have combined the binned results 
for the above three pairs of wavelengths by summing 
the averages for log r 5 4, log r4 2 , log r%\ to get a total 
logr 5 i = log(F A (5100)/F A (1350)). Figure H shows the 
composite color-T max relation for the wavelength range 
5100 - 1350 A as well as this luminosity ratio for individ- 
ual models calculated as described in Sec. 12.31 Over this 
wide wavelength trend, the models show a substantial 
trend toward bluer colors with increasing temperature. 
At low Tmax, there is rough agreement between the ob- 
servations and models, whereas at higher temperatures, 
the observed colors remain fairly constant and depart in- 
creasingly from the models. Richards et al. (2003) find 
a trend of narrower of H/3 with redder continuum slope, 
consistent with our findings. They note the possible con- 
nection with black hole mass (however, they find a trend 
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Fig. 3.— The three continuum flux ratios (5100A/4000A, 
4000A/2200A, 2200A/1350A, as in Figure [2j averaged in bins of 
T max and compared with results of the disk models (assuming 
cos i > 0.5). 
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Fig. 4.— The 5100/1350 A flux ratios for the cos % > 0.5 Kerr 
models together with the composite flux ratios for the SDSS objects 
in T max bins. 

of weaker, not narrower Mg n for the same sample). 

3. DISCUSSION 

We have found that QSOs show mixed agreement 
with the predictions of accretion disk models as a func- 
tion of disk temperature r max . We consider the three 
wavelength pairs rs4, r^, and r2i involving ratios of 
F x at 5100/4000, 4000/2200, and 2200/1350 A, respec- 
tively. Accretion disk models with random inclinations 
show the expected trend of bluer colors with increas- 
ing T m ax- Observed values of r^ agree well with the 
models. For r^, the observations agree with the mod- 
els for Tmax < 10 5 ' 4 K, but for higher temperatures 
the observed colors become redder while the models be- 
come bluer. For ri\ at low T max , the observed colors are 
slightly redder than the models. The observed color re- 
mains roughly constant with increasing T max while the 
models become distinctly bluer. Thus, for both r^ and 
r 2 i, the observed QSOs fail to show the expected trend 
with T max . Possible explanations of this discrepancy 
might include erroneous values of Mbh an d M, host 
galaxy contamination, dust reddening, and inadequacy 
of the disk models. 



3.1. Reliability of disk parameters 

The derived values of T max depend on Mbh and M. 
Values of Mbh from echo measurements of the BLR ra- 
dius (Kaspi et al. 2005, and references therein) agree 
with host galaxy velocity dispersions (Gcbhardt et al. 
2000b; Nelson 2000). Comparison of black hole masses 
from an expression like Eq. |4] with host galaxy luminosi- 
ties and velocity dispersion <r* indicates that Mbh is re- 
liable in the mean to ~ 0.1 dex with a dispersion of 
~ 0.4 dex in individual objects (McLure & Dunlop 2002; 
Greene & Ho 2005). The resulting scatter of 0.2 dex 
and zero-point uncertainty of 0.05 dex in T max does not 
significantly affect our conclusions. 

More difficult to assess is the uncertainty in Al. 
We used our disk models to assess the error intro- 
duced from using a constant bolometric correction /j, = 
Lboi A-La (5100) over our range of T max . The true bolo- 
metric luminosity is given by M and e for the assumed 
black hole spin. We find that the model-derived Jl gen- 
erally exceeds 9 in a way dependent on inclination angle. 
For Kerr holes, taking only inclination angles of less than 
60°, we find values of /l ranging up to ~80. On average, 
our assumed /l of 9 underestimates the true value of /l 
by ~ 0.4 dex. Models with very small cos i can have val- 
ues of Jl of up to ^200-300. Such objects should have 
a much higher T max than we have assigned; but if such 
QSOs can be observed at all, they will represent only a 
small number of objects. Since the bolometric correction 
factor enters into the formula for T max in the one-fourth 
power, the effect of an average 0.4 dex underestimation of 
the luminosity translates into a 0.1 dex underestimation 
of the temperature. 

Figure [5] shows a comparison between the models of 
Fig. [T] and those same models with T max computed from 
Lho\ — 9 x ALa(5100). The effect of underestimating Lboi 
for objects of low inclination is slight, as noted above. 
The effect of significant underestimation of Lboi for high 
inclination objects is evident. These objects' T max are far 
too low compared with their true temperatures. While 
these objects do show a trend of increasing red color with 
T max , this only manifests at low T max and only in objects 
of high inclination. Our discussion has mostly assumed 
a near-extremal Kerr black hole (a* = 0.998). It is not 
known what the black hole spin distribution is in quasars; 
the spin parameter is dependent on the merger history 
of the black holes (which tends to decrease the spin) and 
the gas accretion history (which will spin up the accret- 
ing black hole). It is expected to range from ~0.6 up to 
about a* = 0.9 (see, e.g. Volonteri et al., 2005; Gammie, 
Shapiro, & McKinney, 2004; Hughes & Blandford 2003). 
The effect of assuming a variety of spins for our observa- 
tions would be to lower T max by up 0.46 dex (for a zero 
spin black hole). This may serve to more closely align 
the lower temperature bins with the models, but only 
increases the discrepancy of the observations at higher 
T max . Uncertainties in T max whether due to uncertain- 
ties in Lboi or overestimating the black hole spin do not 
seem likely to resolve the discrepancy between observed 
colors and disk models. 

3.2. Host galaxy contamination 

Another contributor to QSO colors is starlight from 
the host galaxy. The contribution of host-galaxy light 
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Fig. 5. — The same models as computed for Fig. [T]but with T m ax 
computed from the model spectrum using Lbol = 9 x AL^(5100). 
Disks with high inclination show significant under-estimation of 
T max (see discussion of in text). Low inclination disks, however, 
retain the qualitative color v. T max trend seen when the model 
colors are plotted as a function of the true T max . 

is indeed seen via stellar absorption lines in SDSS com- 
posite spectra (Vanden Berk et al. 2001). They find 
a contribution to the composite quasar continuum from 
stars of about 7%-15% at the locations of Ca II A3933 
and Na I A5896 and about 30% at the locations of Ca II 
AA8498, 8542. Studies of black holes in nearby quiescent 
galaxies find Mbh roughly proportional to galaxy lumi- 
nosity (see review by Kormendy & Gebhardt 2000). If 
we assume L ga i oc A/bh, we expect an increasing galaxy 
contribution with decreasing L/L-^d- We show in Fig. [7] 
that for our SDSS QSOs, L/L-Ed increases with T max , so 
the galaxy contribution should make the observed QSOs 
increasingly redder than the models for decreasing T max , 
opposite to the trend that we find. Finally, we have made 
composite spectra for our QSO sample grouped accord- 
ing to T max . The Ca K line, visible but weak, suggests 
a galaxy contribution of only a few percent at A4000, 
and no perceptible difference between the objects with 
higher and lower T max . This suggests that host galaxy 
starlight does not cause the divergence of observed and 
model colors with increasing T max . 

3.3. Reddening 

Numerous studies have considered the degree of red- 
dening of AGN and its wavelength dependence (e.g., 
Wilkes et al. 1999; Richards et al. 2003). Hopkins 
et al. (2004) studied the reddening in QSOs of SDSS 
DR1. Based on correlations of color and spectral curva- 
ture and reddening curves by Pei (1992), they conclude 
that the minority of objects that do show significant red- 
dening are consistent with an SMC reddening curve but 
not with LMC or Galactic reddening. However, the pat- 
terns in Fig. [3] are not consistent with reddening by dust. 
The observed values of 7*42 follow the Kerr models below 
Tmax = 10 K, then turn up sharply toward redder val- 
ues. In contrast, r^x deviates from the models even at the 
lowest temperatures, showing roughly constant values as 
a function of T max . For SMC dust, the reddening of the 
flux ratio is relatively small and inconclusive for our 
discussion. The argument is even stronger for LMC or 
MWG dust, for which the reddening between A1350 and 
A2200 is smaller or even negative. The Gaskell et al. 



(2004) reddening curve is flat in the ultraviolet, giving 
similar values of differential extinction of 7-21 and r^, 
and much less for ^4. Such a reddening curve cannot 
explain the results in Fig. [3l 

We caution that the amounts of reddening required to 
affect the curves in Fig. [3] is not large. An increase in r^i 
of - 0.09 dex requires E{B - V) = 0.03 or A v = 0.09; 
the same change in 7-42 requires E(B — V) = 0.05 or Ay — 
0.16. For comparison, the amount of Galactic reddening 
of the SDSS QSOs used here averages about A g = 0.13 
for the SDSS g-band photometric color. However, the 
reddening would need to correlate systematically with 
broad line width to explain the trends in r42 or 7*21. To 
test this possibility, we assume that the amount of dust 
reddening correlates with broad line width (i.e. T max ) in 
such a way as to give overall agreement in the ratio r%x. 
Figure [B] shows the binned points of Fig. [3] de-reddened 
with an assumed reddening of 7*21 given by A log 7*2 1 = 
—0.57 + 0.125 • logT max . This is tailored in such a way 
as to give overall agreement between the de-reddened 
observations and the models for 7*21,. The SMC extinction 
curve (Pei 1992) gives A log 7-42 /A log r 2 i = 0.61 and 
A log r 5 4/A log r 2 i = 0.16. The figure shows that the 
de-reddened r42 now is bluer than the models for lower 
Tmax, and at high T max it rises to redder values in a way 
not reflected in the models. The de-reddened ^4 now 
is significantly bluer than the models for most values of 
Tmax- Evidently, reddening cannot reconcile the models 
and observation in detail, even for a carefully tailored 
increase of reddening with broad line width. 

Reddening may nevertheless have a significant effect on 
the observed colors and their relationship to the model 
colors. For example, Constantin & Shields (2005) find 
evidence for greater reddening in Narrow Line Seyfert Is 
(NLSls) than in the AGN population as a whole. Baskin 
& Laor (2005) plot the C iv/H/3 flux ratio against con- 
tinuum slope and find evidence for reddening amounting 
to several tenths dex in the line ratio. Even if it does 
not correlate significantly with T mayi , reddening of this 
magnitude would significantly alter the relationship of 
the observed colors and disk models, as in our Fig. [3l 
Willott (2005) argues that the reddening curve proposed 
by Gaskell et al. is affected by selection biases relating 
to the composite spectra used. Higher redshift objects 
define the short wavelength part of the composite, and 
there may be a selection for lower extinction at higher 
redshift because of the limiting magnitude of the survey 
providing objects for the composite. Such a bias might 
affect our results, which use higher redshift objects for 
T2\. However, this bias would cause the extinction to be 
lower for r^i than 7-42, opposite to the result in Fig.[3l 

Baker (1997) finds that the radio-loud quasars (RLQs) 
are redder for more edge-on inclinations. Wills & Browne 
(1986) find broader emission lines in RLQs seen more 
edge-on, as would occur if H/3 lines were emitted from 
a disk-like structure. This would result in broader lines 
for more edge-on disks, and greater reddening for broader 
lines. This is opposite to the trends in Fig. [3] While most 
SDSS QSOs are radio quiet, the same reasoning should 
apply. The relationship between nuclear gas supply and 
luminosity is uncertain. Plausibly, more abundant gas 
might lead to greater reddening as well as greater fueling 
and higher values of Lboi and T max . 
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Fig. 6. — The binned observed points of Fig. [3] de-reddened 
with an SMC extinction curve. We assume here that the difference 
in the 2200 A/f350 A flux ratio between the observed points and 
the models is accounted for by reddening correlated with broad 
line-width. The other two flux ratios are then de-reddened by the 
corresponding amount according to the SMC extinction curve (see 
text) . 
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Fig. 7. — The Eddington ratio (in log). For higher redshift 
objects where A5100 is inaccessible, we scale the 4000 A continuum 
to 5100 A as a power law with a v = —0.44 ± 0.1 (Vanden Berk et 
al. 2001). The uncertainty has only a ~ 1% affect on the derived 
M BH and L hol /L Ed . 



3.4. Disk physics 

Does accretion disk physics underly the difference be- 
tween observed colors and the models? Figure [7] shows 
the strong correlation of T max and iboi/iEd for our sam- 
ple, largely reflecting the fact that lower T max objects 
have broader lines, and hence larger black holes. (Note 
that the derived £boi/^Ed is proportional to the bolo- 
metric correction Jl that we have assumed). Collin & 
Hure (2001) and Collin et al. (2002) have discussed 
high Eddington ratios in AGN based on objects with 
measured reverberation masses for the black hole. They 
note the possible role of slim disks or non-standard disks 
in super-Eddington accretion scenarios. The largest ac- 
cretion rates correspond to the narrowest broad lines, 
particularly in NLS1 objects. In Fig. [31 the common 
NLS1 definition FWHM < 2000 km s" 1 corresponds to 
logT max >5.6. 

Comparing Fig. [7] with Fig.O we see that the deviation 



of from the models occurs for -Lboi/^Ed > 0.3 absent 
significant reddening (see above). This is just the range 
where the inner disk is severely thickened by radiation 
pressure, and departures from a thin disk geometry and 
atmospheric structure might be expected. The concept 
of "slim" accretion disks (Abramowicz et al. 1988) may 
therefore be relevant. These models modify the standard 
equations for thin accretion disks (Shakura & Sunyaev 
1973) by including terms allowing for the radial pres- 
sure gradient and radial advection of heat. Szuszkiewicz, 
Malkan, & Abramowicz (1996) discuss the application of 
slim disks to AGN energy distributions, including the 
soft X-ray excess often observed in AGN (e.g., Turner 
&: Pounds 1989). Slim disks are thin at large radius, 
and then increasingly depart from the thin disk solution 
at smaller radii as the local radiation flux counteracts 
the vertical gravity and thickens the disk. Slim disks 
thicken at a transition radius r s ii m that increases with 
iboi/iEd- Thus, longer wavelength radiation, emitted 
from the cooler, outer disk may obey the thin disk so- 
lution, while shorter wavelengths originating at smaller 
radii depart from the thin disk prediction. We attempted 
to estimate the effect of a slim accretion disk by apply- 
ing an inner cutoff to the disk at the radius where the 
ratio of disk height to radius is H/R « 0.5 (see discus- 
sion of H/R in Hubeny et al. 2000 and Laor & Netzer 
1989). This is based on the idea that the locally ra- 
diated flux and effective temperature may be depressed 
at radii where advection carries energy inward. The L\ 
(5100)/La (4000) ratios are relatively unaffected by an 
inner disk cutoff; however, the color ratios, and r 2 i 
do show an upward trend at higher T max . Although they 
do not reproduce the correct magnitude of the ratios, 
and r 2 i show a flat or 'U'-shaped curve. The result 
of this simple approximation suggests that an accurate 
simulation of slim accretion disks could illuminate our 
observations of objects with high Eddington ratio. 

Comptonization can affect the observed spectra of 
AGN, involving either an hot corona overlying the disk or 
electron scattering in the atmosphere itself (Shang et al. 
2005; Hubeny et al. 2001: and references therein). The 
disk models here do not include Comptonization. Blaes 
et al. (2001), in a study of accretion disk models of 3C 
273, find that Comptonization in the atmosphere does 
not affect the optical and mid-UV spectrum. Neverthe- 
less, Comptonization should be kept in mind in studies 
of the QSO energy distribution. 

We have considered only the optical and ultraviolet 
energy distribution. The statistical sampling allowed by 
SDSS makes it possible to study fairly subtle trends in 
the colors with T max . Nevertheless, T max has a much 
larger effect on the flux at short wavelengths. Laor et al. 
(1997) discuss the soft x-ray properties of AGN. They 
find that accretion disk models, for reasonable distribu- 
tions of -Mbh and M, predict a larger range of optical to 
soft x-ray slope a ox than is observed. Our Figures [3] and 
[4] likewise show a smaller range of color in the observa- 
tions than in the models. 

Boroson and Green (1992) identified a number of cor- 
relations between observable properties of QSOs, the 
strongest of which they designated "Eigenvector 1" . This 
includes a tendency for stronger permitted Fe II broad 
line emission together with weaker narrow [O in] lines. 
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Boroson and Green suggested that Eigenvector 1 was 
driven by Eddington ratio, in the sense of increasing 
Fe II with increasing -Zjboi/^Ed- This relationship was 
confirmed by Boroson (2002), but the physical causes re- 
main unclear. The correlation of iboi/iEd with T max 
suggests that T max may play a role in the physics under- 
lying Eigenvector 1. 
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